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Abstract 

The standing quasi-modes of the ideal MHD in a zero-/? cylindrical magnetic flux tube 
that undergoes a longitudinal density stratification and radial density structuring is consid- 
ered. The radial structuring is assumed to be a linearly varying density profile. Using the 
relevant connection formulae of the resonant absorption, the dispersion relation for the fast 
MHD body waves is derived and solved numerically to obtain both the frequencies and damp- 
ing rates of the fundamental and first-overtone, k — 1,2 modes of both the kink (m = 1), and 
fluting (m = 2) waves. Where k and m are the longitudinal and azimuthal mode numbers, 
respectively. 
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1 Introduction 



Transverse oscillations of coronal loops were first identified by Aschwanden et al. (f999) and 
Nakariakov et al. (1999) using the observations of TRACE. Nakariakov et al. (1999) reported 
the detection of spatial oscillations in five coronal loops with periods ranging from 258 to 320 s. 
The decay time was (14.5 ± 2.7) minutes for an oscillation of (3.9 ± 0.13) millihertz. Also Wang 
& Solanki (2004) described a loop oscillation observed on 17 April 2002 by TRACE in 195A. 
They interpreted the observed loop motion as a vertical oscillation, with a period of 3.9 minutes 
and a decay time of 11.9 minutes. All these observations indicate strong dissipation of the wave 
energy that may be the cause of coronal heating. 

Since the discovery of the hot solar corona about 66 years ago, different theories of coronal 
heating have been put forward and debated. Ionson (1978) was first to suggest that the resonant 
absorption of MHD waves in coronal plasmas could be a primary mechanism in coronal heating. 
Since then, much analytical and numerical work has been done on the subject. Rae and Roberts 
(1982) investigated both eikonal and differential equation approaches for the propagation of 
MHD waves in inhomogeneous plasmas. Hollweg (1987a,b) considered a dissipative layer of 
planar geometry to study the resonant absorption of coronal loops. Poedts et al. (1989, 1990) 
developed a finite element code to elaborate on the resonant absorption of Alfven waves in 
circular cylinders. 

Davila (1987) and Steinolfson & Davila (1993) did much work on resonant absorption through 
resistivity. Goossens et al. (2002) used the TRACE data of Ofman & Aschwanden (2002) to 
infer the width of the inhomogeneous layer for 11 coronal loops. Ruderman &: Roberts (2002) 
did similar analysis with the data of Nakariakov et al. (1999). Van Doorsselaere et al. (2004a) 
used the LEDA code to study the resistive absorption of the kink modes of cylindrical models. 
They concluded that, when the width of the nonuniform layer was increased, their numerical 
results differed by as much as 25% from those obtained with the analytical approximation. In 
the vicinity of singularity, field gradients are large. Recognizing this, Sakurai et al. (1991a, b) 
and Goossens et al. (1992, 1995) developed a method to analyze dissipative processes in such 
regimes and to neglect them elsewhere. 

Safari et al. (2006) studied the resonant absorption of MHD waves in magnetized flux 
tubes with a radial density inhomogeneity. Using the approximation that resistive and viscous 
processes are operative in thin layers surrounding the singularities of the MHD equations, they 
obtained the full spectrum of the eigenfrequencies and damping rates of the quasi-modes of the 
ideal MHD of the tube. Both surface and body modes were analyzed. 

Verwichte et al. (2004), using the observations of TRACE, detected the multimode oscil- 
lations for the first time. They found that two loops are oscillating in both the fundamental 
and the first-overtone standing fast kink modes. According to the theory of MHD waves, for 
uniform loops the ratio of the period of the fundamental to the period of the first overtone is 
exactly 2, but the ratios found by Verwichte et al. (2004) are 1.81 and 1.64. However, these 
values were corrected to 1.82 and 1.58, respectively, by Van Doorsselaere et al. (2007) and thus 
clearly differ from 2. This may be caused by different factors such as the effects of curvature 
(see e.g. Van Doorsselaere et al. 2004b), leakage (see De Pontieu et al. 2001), magnetic twist 
(see e.g. Erdelyi &: Fedun 2006; Erdelyi &; Carter 2006) and density stratification in the loops 
(see Andries et al. 2005b; Erdelyi & Verth 2007). However, now it clear that the only cause 
of this deviation is the longitudinal density stratification. Curvature does not affect sufficiently 
the frequencies. Leakage can cause the damping of oscillations but practically does not affect 
the frequencies. Recently Ruderman (2007) showed that the magnetic twist also does not affect 
the frequencies. 
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Andries et al. (2005b) considered a coronal loop model with a straight cylindrically magne- 
tized flux tube. They elaborated the effect of longitudinal density stratification on the oscillation 
frequencies and the damping rates of fast surface waves by resonant absorption. 

Andries et al. (2005a), with the help of the method introduced by Andries et al. (2005b), 
showed that in the presence of a longitudinal density stratification the ratio of the periods of 
the fundamental kink mode of a coronal loop and of its first harmonic is lower than 2. They 
used this ratio to estimate the density scaleheight in the solar atmosphere. 

Arregui et al. (2005) studied the effects of both radial and longitudinal density stratifications 
on the frequency and damping rate of resonantly fast kink modes. Using a numerical code 
(POLLUX), they solved the linear resistive MHD equations for a cylindrical flux tube model 
with a wide range of values of several loop parameters. 

Safari et al. (2007) investigated the oscillations of coronal loops in presence of vertical 
stratification for a zero-/? plasma. They solved the radial equation in the thin tube approximation 
and the transverse equation by both perturbational and numerical techniques. They confirmed 
the result obtained by Andries et al. (2005a) that the ratio of periods of the fundamental and 
first-overtone modes differ from 2 in loops with longitudinal density stratification. 

Karami & Asvar (2007, hereafter Paper I) studied both the oscillations and damping of 
standing fast MHD body waves in coronal loops in the presence of both longitudinal density 
stratification and radial density structuring. The radial structuring was assumed to have a step- 
like density profile. They derived both the frequencies and damping rates of the fundamental, 
first-overtone and second-overtone frequencies of both the kink and fluting modes. They obtained 
that the ratios of the frequencies of the first-overtone and its fundamental mode for both the 
kink and fluting modes are lower than 2 (for unstratified loops). 

In the present work, our aim is to investigate the effects of longitudinal density stratification 
and radial density structuring on resonant absorption of standing quasi-modes of the ideal MHD 
in the cold coronal loops observed by Verwichte et al. (2004) deduced from the TRACE data. 
This paper is organized as follows. In Sections 2 and 3 we combine the two techniques of Paper 
I and of Thompson &: Wright (1993) to derive the equations of motion, introduce the relevant 
connection formulae and obtain the dispersion relation. In Section 4 we give numerical results. 
Section 5 is devoted to conclusions. 



2 Equations of motion 

The linearized MHD equations for a zero-/3 plasma are 

^ = -^{(V x <5B) x B + (V x B) x SB} + ^V 2 <5v, (1) 
ot 4irp p 

dSB _ ,„ . c 2 2 



« =Vx(*vxB) + — V'ffl, (2) 

where <5v and 5B are the Eulerian perturbations in the velocity and magnetic fields; p, a, rj 
and c are the mass density, the electrical conductivity, the viscosity and the speed of light, 
respectively. The simplifying assumptions are the same as Paper I. 

From Andries et al. (2005b), in the absence of dissipations, the perturbed quantities 5~B(r, z) 
and 5v(r, z) can be expressed in terms of a complete set of orthonormal functions ip^(z) as 
follows 

oo 

SB{r,z) = Y / SB^(r)^ k \z), (3) 
k=i 
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where the same relation is kept for <5v(r, z). Note that ijj( k \z) itself satisfies the eigenvalue 
relation L,Aip( k \z) = \kip^ k \z), where La is the Alfven operator, 

2 B 2 d 2 / 2 2 # x 

with Alfven velocity va = an d straight constant background magnetic filed B = Bz. 

One may show that the resonance absorption occurs where the Alfven operator has a vanishing 
eigenvalue. 

Let us denote the radius of the tube by R and a radius at with the resonant absorption 
occurs by R\ < R. The thickness of the inhomogeneous layer, I = R — R\, will be assumed to 
be small. From Safari et al. (2006) and Paper I, the density profile is assumed to be 



71" Z 

p(r,z) = p (r)exp [ - ^sin( — ) 



Po(r) = < 



R-Rl 



Pin, (r < Ri), 

{R-r) + p cx , (R 1<r <R), (5) 
Pex, {r>R), 



where p is defined as stratification parameter, H and L are the density scale hight and length 
of the loop, respectively. Also p- m and p cx are the interior and exterior constant densities of the 
tube, respectively. 

From paper I, in the absence of dissipations, in the interior region, r < R\, solutions of Eqs. 
(l)-(2) are 

{I m {\hn,k\r), k 2 n k<0, surface waves, 
J m {\hn,k\r), fci 2 n , fe >0, body waves, (6) 
U2 _ V.fc 
^in.fe — B 2 /4vr ' 

where J m and I m are Bessel and modified Bessel functions of the first kind, respectively. In the 
exterior region, r > R, the waves should be evanescent. The solutions are 

«?(«.*>(r) = K m (k cx , k r), k 2 e ^ k = -^jf- > 0, (7) 

where K m is the modified Bessel function of the second kind. 

Here we only consider the body waves. Because under coronal conditions, VA in < ^A ox5 
magnetic flux tube supports fast body oscillations and there are no longer any surface modes 
(see Fig. 4 in Edwin k, Roberts 1983). From Eq. (7) in Paper I and Eqs. (3), (6), one can 
obtain 5Bz n ' k \r, z) and 5vr U ' k \r, z) in the interior region (r < R±) as 

+oo 

6B<r>(r,z) = Y. AMj rn{\hn,k\r)^ k) {z), 



k=l 



6vW{r,z) = -^E^^ (m ' fe) ^(l^|r)^(,), 



(8) 



where a prime on J m and hereafter on each function indicates a derivative with respect to their 
appropriate arguments. The results for the exterior region are the same as Eq. (8), except that 
J m and index "in" are replaced by K m and "ex" everywhere. 
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3 Connection formulae, dispersion relation and damping 



From Andries et al. (2005b) and according to the connection formulae developed by Thompson 
& Wright (1993), the jump across the boundary for 5B Z and 5v r are 

[6B Z ] = 0, (9) 



[SVr] = "E 



5B 



(in,fc) 



k=i ir\^ k )\L Al 



(in,fe) 



0(in,fc) j 



where 



L A (/) iin ' k) = = 



^(in,fc) 



and 



dL A 



dr 



( PinCl + SfcOU 2 -^) ^ J 
1 1 k=j 

dp {r) 



= UJ (1 + Sfcfc, 



r=r A 



2 f L kir r . tt i . jvi 



(10) 

(11) 
(12) 

(13) 
(14) 



R — R\ 



-L0 2 (l + S kk ). 



(15) 

Note that R\ < r A < R is the radius at which the singularity occurs and, Co = co + i"f where 7 
is damping rate. Substituting the fields of Eq. (8) in jump conditions, gives 



/ n' 



jj(ex,l) _jjl m ' 1 ^ jjl ex ^> _-\j( ln ' z > 

~(ex,l) „(in,l) ^(in.1) „(ex,2) „(in,2) ^( 



(in,l) 



(ex,2) 



(in,2) 



in,2) 



|-r(ex,l) 

„(ex,l) 
^2 



-n 

j(in,l) 
J 2 



(in,l) 



n 



(ex,2) 



n (in,2) 



p (in,l) s (ex,2) ^(in.2) + p(in,2) 



^(in,l) 
^(ex,2) 
^(in,2) 



0, 



(16) 



/ 



where 



n 



(in,fc) 
3 

(ex,fc) 



n 

j 

^(in.fe) 
"J 

— ,(ex,fc) 



p(in,fc) 



JjJ^kj ,(in,fc) 
g^fc) ^(ex,fc) 



2^,2 A( in ' fe )rj( in ' fe ) 



z=i 



,(in,fc) 



4R3/^(in,fc) 



(in,fc) 



(17) 
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and 



(in,fc) 



J \ttR> j 



2^? 



k = j 



(18) 



4^ = { (^)(^) 2 ^ ^i, (19) 

1 k — j 

Xk = \km,k \ R an d Uk = k cx ,kR- Note that the dispersion relation is then given by requiring that 
the system (16) has non-trivial solutions i.e. its determinant is zero. 



4 Numerical results 

As typical parameters for a coronal loop, we adopt radius = 10 3 km, length = 10 5 km, p m = 
2 x 1CT 14 gr cm~ 3 , p C x/Pin = 0.1, B = 100 G. For these parameters one finds v Ain = 2000 km 
s" 1 , v Ac ^ = 6400 km s" 1 and uj A . n := ^ = 0.02 rad s -1 . 

The effects of density stratification on both the frequencies co and damping rates 7 are 
calculated by numerical solution of the dispersion relation, Eq. (16). The results are displayed 
in Figs. 1 to 7. Figures 1-4 show the frequencies and damping rates as well as the ratio of the 
oscillation frequency to the damping rate of the fundamental and first-overtone, k = 1,2, kink 
(m = 1) and fluting (m = 2) modes versus the stratification parameter p = L/ttH for a loop 
with constant length and varying scale hight. Figures 1 to 4 reveal that: i) Both frequencies, 
0J1, 0J2 and their corresponding damping rates, |7i|, [72) increase when p increases. The results 
for uj\ and 71 are in agreement with those obtained by Andries et al. (2005b) and Arregui et 
al. (2005). Also the results for all uis and 7s are in accordance with that obtained in Paper I. 
ii) For m = 1 with increasing p, the damping rates, for instance I71I, increase (~ 2-170 percent) 
compared with a non-stratified loop. 3) The ratio of the oscillation frequency to the damping 
rate, w/l7l> ^ s independent of stratification. This result is in good concord with those obtained 
numerically by Andries et al. (2005b) and analytically by Dymova & Ruderman (2006). 

Figure 5 same as Fig. 1 shows the result of frequency, damping rate and ratio uj/\j\ for 
the fundamental kink modes but for a thick inhomogeneous layer with l/R = 0.1. Comparing 
the results obtained for the ratio W1/I71I in Fig. 1 and 5 shows that the values of uj/\j\ for 
l/R = 0.02 in Fig. 1, five times greater then the case of l/R = 0.1 in Fig. 5. This result is in 
good agreement with that obtained by Dymova & Ruderman (2006) who showed that under the 
assumption of uniform stratification, i.e. p cx (z) / pj n (z)=const., the ratio u)/\j\ for kink modes in 
thin tube-thin inhomogeneous layer approximation is proportional to (l/R) -1 and is independent 
of stratification. 

The ratio of the frequencies 002/^1 of the first-overtone and its fundamental mode for both the 
kink (m = 1) and fluting (m = 2) modes are plotted in Fig. 6. Figure shows that for both modes, 
the ratio of the frequencies decreases from 2 (for unstratified loop) and approaches below 1.4 
with increasing density stratification. The result of U2/0J1 for m = 1 is in accordance with that 
obtained by Andries et al. (2005a) and Safari et al. (2007). Also the results for both m = 1 and 
m = 2 are in agreement with those obtained in Paper I. For m = 1, p = 0.59 and 0.89, the ratios 
U2/0J1 are 1.821 and 1.576, respectively. These are in good agreement with the frequency ratios 
observed by Van Doorsselaere et al. (2007), 1.82±0.08 and 1.58±0.06, respectively, deduced 
from the observations of TRACE. Fig. 6 does not show a noticeable difference between the kink 
and fluting modes, but our numerical values shows that for a given stratification parameter p, 
the frequency ratio U2/0J1 increase slightly when the azimuthal mode number m increases. 
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Figure 7 shows the ratios of frequencies co^^/ujf 1111 and damping rates ^thick^thm Q f 
fundamental, k = 1, kink modes (to = 1) versus the stratification parameter for two tube radii 
R/L = 0.02 (thick) and 0.01 (thin). Figure 7 clears that both the frequencies and damping rates 
are independent of the tube radius, R, in the limit of slender tubes (see e.g. Van Doorsselaere 
et al. 2004a). 

5 Conclusions 

Resonant absorption of standing fast MHD body waves in coronal loops in presence of both 
the longitudinal density stratification and radial density structuring is studied. The radial 
structuring is assumed to be a linearly varying density profile. To do this, a typical coronal loop 
is considered as a straight pressureless cylindrical flux tube embedded in a constant background 
magnetic field. Using the relevant connection formulae, the dispersion relation is obtained and 
solved numerically for obtaining both the frequencies and damping rates of the fundamental and 
first-overtone kink and fluting modes. Our numerical results show that: 

i) Both frequencies and damping rates of the fundamental and first-overtone modes of both 
the kink (to = 1) and fluting (to = 2) waves increase when the stratification parameter increases. 

ii) The ratio of the oscillation frequency to the damping rate, for both the kink (to = 1) 
and fluting (to = 2) modes is independent of stratification. 

iii) The ratio of the frequencies u^A^i for both the kink (to = 1) and fluting (to = 2) modes 
are lower than 2 (for unstratified loops), respectively, in presence of the longitudinal density 
stratification. The result of w^/oji for kink modes is in accord with the TRACE observations. 

Note that in our work, calculating the frequencies 0J\,0J2 and their corresponding damping 
rates 71,72 as well as the ratio of frequencies u^/Vl, for fluting modes (m=2) is a new result in 
comparison with Andries et al. (2005a,b) and the rest is a repetition of the results previously 
obtained by other authors. 
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Kink Modes (m=l) 




^ 2 
o 



— 1.5- 

: 

' 1 ' 1 1 1 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
H=L/tiH 

Figure 1: Frequency of the fundamental kink mode (m = 1) and its damping rate as well as the 
ratio of the oscillation frequency to the damping rate as a function of the stratification parameter 
fi = L/ttH, for a loop with constant length and varying scaleheight. The loop parameters are: 
L = 10 5 km, R/L = 0.01, l/R = 0.02, p e / Pi = 0.1, Pi = 2 x 10~ 14 gr cm' 3 , B = 100 G. Both 
frequencies and damping rates are in units of the interior Alfven frequency, uja { = 0.02 rad s -1 . 
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Kink Modes (m=l) 




H=L/tiH 

Figure 2: Same as Fig. 1, for the first-overtone kink modes. 



Fluting Modes (m=2) 
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x. 0.8 



3 0.6- 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
H=L/tiH 

Figure 3: Frequency of the fundamental fluting mode (m = 2) and its damping rate as well 
as ratio of the frequency and the damping rate as a function of the stratification parameter 
H = L/ttH. Auxiliary parameters as in Fig. 1. 
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Fluting Modes (m=2) 
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Figure 4: Same as Fig. 3, for the first-overtone fluting modes. 



Kink Modes (m=l) 




H=L/tiH 

Figure 5: Same as Fig. 1, for the fundamental kink mode (m = 1) with l/R = 0.1. 
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Kink Modes (m=l) 




Figure 6: Ratio of the frequencies 002/^1 of the first-overtone and its fundamental mode versus 
/U = L/ttH for both kink (to = 1) and fluting (m = 2) modes. Auxiliary parameters as in Fig. 
1. 
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Figure 7: Ratios of the frequencies uf 11 ^ / uf 1111 and the damping rates ^ hlck /^ hm of the funda- 
mental kink (to = 1) mode versus fi = L/ttH for two tube radii R/L = 0.02 (thick) and 0.01 
(thin). Other auxiliary parameters as in Fig. 1. 
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